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Abstract 

Models with large extra dimensions as well as unparticle models could give rise to new 
phenomena at collider experiments due to real emission or virtual exchange of gravitons 
or unparticles. In this paper we present the common implementation of these processes in 
the Monte Carlo generator Pythia8, using relations between the parameters of the two 
models. The program offers several options related to the treatment of the UV region of the 
effective theories, including the possibility of using a form factor for the running gravitational 
coupling. Characteristic results obtained with Pythia8 have been used to validate the 
implementations as well as to illustrate the key features and effects of the model parameters. 
The results presented in this paper are focused on mono-jet, di-photon and di-lepton final 
states at the LHC. 



1 Introduction 



Models with large extra dimensions (LED) are popular extensions of the Standard Model (SM) 
and studies of LED phenomena are usually based on the so-called ADD scenario [1], where 
gravity alone has access to the extra dimensions. The large size of the extra dimensions gives rise 
to dense Kaluza-Klein (KK) mass modes which appear as a continuous graviton mass spectrum. 
Due to the large number of contributing mass states graviton processes could obtain sufficiently 
large cross sections to be observed at collider experiments, where the new processes involve real 
emission or virtual exchange of gravitons. These processes have been studied in great detail [2], 
however, no signs of new physics have been indicated at existing experiments. 

More recently so-called unparticle models [3] have also gained a large interest. These mod- 
els relate to physics originating from a scale invariant new sector which is coupled to the SM 
through a connector sector with a high mass scale. Unparticle models can often, at least from a 
phenomenological point of view, be associated with extra dimension models and the particular 
case addressed here would result in phenomena very similar to the scenario of LED gravity. 
Having access to processes with the same final state from more than one model can be helpful to 
generalise the experimental analysis as well as for studying differences in the measured observ- 
ables. The latter, in turn, could also be of value after a discovery to characterise the observed 
signal. 

This paper describes the processes involving real emission or virtual exchange of unparticles 
(U) or LED gravitons (G) as implemented in Pythia8 [4]. The corresponding processes are: 

• Mono-jets (U/G + jet): gg ->U/G + g, qg-^U/G + q, qq-^U/G + g; 
. Mono-Z/7 (U/G + Z/j): ff^U/G + Z, ff^U/G + T , 

• Di-photons (77): // ->■ U*/G* -)• 77, gg -)• U*/G* ->■ 77; 

• Di-leptons (££): ff -»• U*/G* -> ££, gg -> U* /G* -»• ft. 

Because of the large similarities both models can be covered by a common implementation [5], 
where only a simple translation of a few model related constants is required to change between 
the models. In addition to technical advantages, this common description makes comparisons 
between the processes both consistent and transparent. The effective theory description of these 
processes behaves poorly in the UV limit and several related options are available in the program. 
These include a form factor [6] for the graviton processes which represents a realistic alternative 
of the behaviour at high energies due to the running of the gravitational coupling. 

The Pythia8 implementations have been validated against the original theory papers and we 
present here some of these results. Besides providing reference results for the Pythia8 program, 
the results illustrate different characteristic effects of the model parameters and various key 
features 1 . All results are produced with respect to the nominal LHC conditions (pp-collisions at 
a centre-of-mass energy of y/s = 14 TeV), unless explicitly stated otherwise. 

This paper is organised as follows. The first section starts with a description of the model 
parameters and different conventions. It is followed by a description of the implementation 
of the processes and the relations between the two models. This section also describes the 
available options for treating the UV region of the effective theories. Sec. 3 to 5 present the 

1 In order to prevent clutter, the figures presented in this paper do not include statistical errors, which should 
not be significant for the conclusions. 
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different program options in Pythia8 together with results related to the individual processes: 
14 /G + jet, 77, and U production. In addition, Sec. 5 draws special attention to possibile effects 
on the forward-backward asymmetry in lepton pair production from spin-1 unparticle exchange. 
The U/G + Z/j processes are not addressed here since they were already presented in a previous 
paper [5]. A complete list of the process and model parameter names used in the Pythia8 
program can be found in the appendix. The appendix also contains complementary information 
used for the scalar unparticle process, // — > Z14, which is not available in the literature. 

This paper presents results obtained with Pythia version 8.125 and all parameters not 
explicitly mentioned in the text have been kept at the default values of this version. During the 
work a few minor corrections have been made 2 , which all are included in the following version 
8.130. In order to produce results comparable to what is available in the literature, only the 
hard part of the process has been considered 3 and the MRST2001 parton distribution functions 
(PDF) of the proton have been used unless stated otherwise. 

2 Implementation in Pythia8 

The common implementation of the LED graviton and unparticle processes in Pythia8 is gen- 
erally based on the cross section expressions derived for the unparticle case. These expressions 
can be converted into the corresponding LED graviton formulae by simple translations of a few 
model related constants. For this reason we first discuss the unparticle process and afterwards 
we introduce the translation into its LED equivalent. The values of the model parameters are 
sometimes constrained by arguments valid in different scenarios of the model, however, in order 
to allow for more phenomenologically motivated studies the implementations in PYTHIA8 gen- 
erally allow for any parameter value that is technically possible. Particular model constraints 
therefore have to be respected explicitly by the users. As mentioned above, the common imple- 
mentation of the U/G + Z/j processes was originally presented in [5], together with a study of 
U/G + Z production at the LHC, but for completeness relevant parts related to these processes 
will be repeated here. 

2.1 Parameter Conventions 

The parameter convention used here for the unparticle case follows [7]. It is assumed that any 
potential gaps in the unparticle mass spectrum are sufficiently small not to affect the high energy 
behaviour at particle colliders and that the life time is long enough to prevent decays at distances 
relevant for collider experiments. The main parameters are a scale dimension parameter {du), 
an unparticle renormalisation scale (A^) and a coupling constant (A) of the unparticles to the 
SM fields, which is related to the mass scale of the connector sector. 

In principle there are restrictions on the allowed ranges of the parameters. Full conformal 
invariance restricts the scaling dimension of the unparticle to du > 1 (scalar), du > 3 (vector) 
and du > 4 (tensor), which follow from Mack's unitary constraint, as discussed for example 
in [8]. By assuming scale invariance, but not full conformal invariance, these restrictions are 
relaxed. In this paper we take a phenomenological approach and allow for the full du range and 

2 The corrections are related to the 11 process and, in addition, the spin-0 case of the U + Z/j process has been 
added. 

3 Where the program parameters PartonLevel:MI, PartonLevel : ISR and PartonLevel :FSR have been turned 

off. 
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the lowest bound du > 1 is used as the overall lower limit in Pythia8. For some processes it is, 
however, necessary to constrain this parameter to the range 1 < dy < 2 for numerical reasons. 
For du > 2, especially in the case of virtual unparticle exchange, contributions that depend on 
the UV completion of the theory are relevant which suggests that the effective theory used will 
not be valid. The ranges allowed in Pythia8 are in accordance with [7, 9], where for example 
Eq. (4) in the latter paper shows that the propagator diverges 4 when du — > 2, 



where Ad u is a normalisation constant related to the rather peculiar unparticle phase space. 
Hence, this divergence implies that the virtual U exchange processes in Pythia8 require du < 2. 
The formulae obtained for the U emission processes on the other hand give finite results also 
for the range du > 2. However, in this case the results have to be interpreted more carefully as 
discussed in Sec. 2.3. 

The parameter convention used here for LED gravity follows [11]. The gaps between the 
individual KK modes of the graviton are determined by the size of the extra dimensions together 
with a possible curvature of space-time, like in the so-called Randall-Sundrum scenarios [12]. 
It has been shown that a small curvature can remove cosmological constraints from low energy 
processes [13]. However, for the processes addressed here it is assumed that any mass gaps in the 
graviton mass spectrum are sufficiently small for the high energy phenomena to be unaffected. 

The main model parameters correspond to the fundamental scale of D (= n + 4) dimensional 
gravity (Mp), the number of extra dimensions (n) and in the case of virtual graviton exchange an 
effective cut-off scale (Ay). The extra dimensions are normally assumed to form a n-dimensional 
torus with a common radius (Rd)- This radius is determined by n, together with the 4- 
dimensional Planck mass and is for this reason not included as a parameter for the processes in 
Pythia8. A number of different conventions are used in the literature and the most popular 
are related as follows: 



• GRW [11]: M D , R D , n, A T (used in Pythia8); 

• MPP [14]: M n+2 = 2M™+ 2 ; 

• HLZ [15]: M£+ 2 = 87r 1_ t T(f )M™ +2 , R s = 2irR D ; 
. HR [6]: Af, = f A 4 . 

2.2 Implementation of the Processes 



The processes in Pythia8 are available for unparticles with spin 0, 1 and 2. However, not all 
spin values are available for all processes. The U/G + Z/*y processes are implemented as 2-to— 2 
parton level processes based on the Z matrix elements (ME) calculated in [7] for spin-1 and 
spin-2 unparticle emission. In addition the spin-0 case is available, which is based on the ME 
given in the appendix. In the case of Z production, no Z/-y* interference is included and the 
Z decays isotropically. The photon process corresponds to the photon limit of the Z process. 

4 The relation between du and the UV sensitivity has for example been studied in [10]. 
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All U/G emission processes implemented in Pythia8 share the graviton particle code 5000039, 
since the U or G has a continuous mass spectrum. This is in contrast to processes related to 
an individual mass state, where for example the first KK mode of the graviton would have the 
particle code 5100039. The main difference with respect to existing processes in Pythia8 was 
the continuous mass spectrum of the unparticle. This has been solved by generating unparticle 
masses, using already existing Pythia8 functionality, according to a Breit-Wigner distribu- 
tion. The events are then re-weighted in order to obtain a correct mass distribution, as given 
by the differential cross section. For this reason the Breit-Wigner shape should be matched 
to the relevant range of the mass spectrum in order to achieve the highest possible Monte 
Carlo (MC) efficiency and generation speed. The Breit-Wigner shape is defined by the program 
parameters, 5000039 :m0, 5000039 :mWidth, 5000039 :mMin, 5000039 :mMax, according to the 
standard particle data scheme in Pythia8. 

The U/G+jet processes have been implemented in the same way as the U/G+Z/j processes. 
They are available for spin-0 and spin-1 unparticles, however, graviton emission is the only 
available spin-2 scenario. The processes are based on MEs from [11] for G emission and [7, 16] 
for the U processes. The spin-0 unparticle interactions used for the processes presented here are 
given by the effective operators [7], 



where / refers to fermions and G a p to gauge bosons. Both terms in Eq. (3) lead to amplitudes 
at leading order for, qq — > U/G + g and qg — > U /G + q, whereas gg — >■ U /G + g only involves 
the gluon vertex, i.e. the second term. Since the gluon vertex is suppressed by one additional 
power of An compared to the fermion vertices, i.e. two powers with respect to the cross section, 
the MEs of the two quark initiated processes are approximated by including only diagrams 
with unparticle emission from a fermion line. The coupling constants of the two terms are not 
necessarily the same and it could be of interest to study them separately. Again due to the 
suppression by different powers of An in the two terms, this can be approximated by using 
qq — > U/G + g and qg — > U/G + q to represent the contributions from the quark vertices and 
gg — > U/G + g to represent the gluon vertex contribution. 

In order to obtain the corresponding graviton emission processes, the unparticle parameters 
have been interpreted as follows, 



where A2 and X' 2 refer to the possibility of having different couplings related to the two effective 
spin-2 unparticle operators [7]. This translation is automatically done inside Pythia8 when the 
graviton processes are used. 

The virtual U/G exchange processes are implemented in P ythia8 as 2-to-2 parton level 
processes. Certain processes contain coherent SM amplitudes and converge to the SM results if 
the unparticle contribution is close to zero, e.g. by setting A = or Ay — > 00. The SM results 





A du -> S'(n) 



r = A' 2 /A 2 = 1 
n 

du = ^ + l 
A u = M D 
, N 2ttW 2 
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shown for comparison in the figures below have been obtained by the corresponding U/G process 
using such extreme parameter values. 

The 77 processes are available for spin-0 and spin-2 unparticles based on the matrix elements 
in [7, 9]. In the spin-2 case, the fermion initiated process includes the corresponding SM t-channel 
process (equivalent to PromptPhoton : f f bar2gammagamma) together with the related interference. 
The spin-2 gluon process does not involve the possible SM box diagram which for this reason 
must be included separately, e.g. using the PYTHIA8 process 5 PromptPhoton :gg2gammagamma. 

The U processes only concern charged lepton pair production. The fermion initiated process 
is available for both spin-1 and spin-2 unparticles and includes the SM Z/j* amplitudes and 
interference. The gluon process is only available for spin-2 and does not contain any SM con- 
tributions. All processes are implemented following the MEs in [7]. Since the fermion initiated 
process could be used for studies at lepton colliders, it should be noted that currently no SM 
t-channel diagram is included in Pythia8. For example, this would be relevant for studies of 
Bhabha scattering. The SM contribution should therefore be similar to the Pythia8 process, 
WeakSingleBoson:ffbar2gmZ, when only the charged lepton final states are considered. A spin- 
1 unparticle can have chiral couplings to fermions which could lead to interesting interference 
effects [17]. For this reason the chiral properties of the coupling can be specified explicitly 
by two coupling parameters (gxx, 9xy)- These can be assigned the values 1, or -1 which 
are multiplied with the different helicity amplitudes, XX = LL/RR and XY = LR/RL. 
This implies that a vector coupling corresponds to, gxx = 9xy = 1, and an axial coupling to, 
9xx = 1,9xy = -i- 

The corresponding graviton processes, where the unparticle is replaced by a virtual graviton, 
are again obtained using the common implementation. This has been achieved by the same 
approach as for the U/G emission processes where two constants had to be adjusted. The cross 
sections given by the spin-2 unparticle MEs discussed above become identical to the graviton 
cross section expressions in [11] when the unparticle parameters are translated as follows, 



The virtual graviton exchange processes are UV sensitive, similar to the unparticle scenario. 
The graviton amplitudes are dominated by contributions that depend on the UV completion of 
the theory when the number of large extra dimensions 6 , n > 2. Since this scenario is generally 
of great interest, the UV sensitive part is commonly replaced by an arbitrary cut-off scale, which 
parametrises our ignorance and removes the n dependence. 

2.3 Treatment of the Effective Theory at High Energies 

Since the current limits on the effective mass scales related to the different processes (A^, Mn 
and At) are well below the possible center-of-mass energies of the hard process, s, at the LHC, 
a number of different options are available for the case where the hard scale of the process 

5 Note that this process is based on a massless approximation which can be relevant. 
6 Which have a similar role to d u in the unparticle model. 



du = 2 
Xu = 1 

A = V^TT 

A u = A T 



(8) 
(9) 
(10) 
(11) 
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approaches or exceeds the scale of validity of the effective theory. These are specified by the 
program parameter CutOf fmode, which can take the following values: 

• CutOffmode = 0: Include all s values; 

• CutOffmode = 1: Suppress the cross section for s > Ay (M|,); 

• CutOffmode = 2: Gravity form factor, using SigmaProcess :renormScale2; 

• CutOffmode = 3: Gravity form factor, using /i = E* et , z ,. 

The default option (CutOffMode = 0) implies no restrictions, but includes the cross section 
contribution for any s value. The first alternative option (CutOffMode = l) simply suppresses 
the cross section [18] by a factor A u /s 2 at s values that exceed the mass scale of the effective 
theory, Ay or M|, . This truncation of the cross section also implies that the U/G mass spectrum 
is suppressed at larger values. The U/G mass spectrum becomes increasingly peaked towards 
larger values with increasing n or du- For this reason, the suppression effect generally becomes 
more significant for increased values of these parameters. In a similar way the effect of this 
truncation becomes larger with an increased transverse energy requirement in the analysis. This 
option is only implemented for the U/G emission processes and further details can be found 
in [5]. Starting the truncation directly above the effective mass scale can often be considered 
as conservative, where for example [11] estimates the effective theory for the LED graviton 
processes to be perturbative up to about 7Mr>. This rather crude method of truncating the 
cross section was primarily intended to be used for verifying that the truncated part, which 
corresponds to the region where the effective theory might not be valid, has a negligible effect 
on the total cross section. However, in the case where it is not negligible the truncated value 
could also be used as a conservative estimate. 

The possibility of using a form factor [6] for the gravitational coupling has also been imple- 
mented, 

-i 

(12) 

Here \x is a renormalisation scale and t is a 0(1) free parameter which relates to the unknown 
details of the running of the gravitational coupling. The parameter M is associated with the 
cut-off scale of the effective theory {Mr, or At depending on the process) . The form factor leads 
to a weaker gravitational coupling at higher energies and this provides a realistic alternative of 
the behaviour when s approaches M 2 . For both real emission and virtual exchange of gravitons 
the renormalisation scale fi can be specified (CutOffMode = 2) to follow the choice made by the 
program parameter, SigmaProcess : renormScale2. For graviton emission it is possible to set 
the renormalisation scale equal to the jet/Z/'j centre-of-mass energy (CutOffMode = 3) which 
was used in [6]. It is argued in [6] that the form factor can generally prevent virtual graviton 
exchange amplitudes from violating unitarity by requiring t < 2, which is suggested as an upper 
bound on this parameter. 

3 Real Emission, Mono-Jets 

Both the LED graviton and unparticle models can lead to final states with a single jet plus 
missing transverse energy, with balancing transverse momenta. The related G emission processes 



F(ji;t,M,n) = 



1 + 
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Figure 1: Differential cross section of G + jet production at the LHC (y/s = 14 TeV) as a function of 
transverse jet momentum. All shapes are normalised to the case n = 2. 



available in PYTHIA8 consist of qqbar2Gg , qg2Gq, gg2Gg, with the program parameters n, MD , 
CutOffmode and t. The corresponding U emission processes are qqbar2Ug, qg2Uq, gg2Ug, 
with the relevant parameters spinU, dU, LambdaU, lambda and CutOffmode. The unparticle 
processes are only available for the spin values spinU = 0,1 and, in both the case of G and U 
emission, the user should match the generated mass distribution to the cross section, as discussed 
in section 2.2. 

This section presents characteristic results related to U /G + jet production at the LHC 
obtained using Pythia8. The processes are validated with respect to the results presented in 
[7, 11, 16] and, to allow for easy comparison, most of the selection cuts and model parameter 
values have been used consistently with these papers. 



3.1 LED Graviton Emission 

In Fig. 1 the differential cross section of G + jet production at the LHC is shown as a function 
of the transverse jet momentum for different values of Md and n. For these generator level 
results the jet simply refers to the outgoing quark or gluon in the hard interaction. The jet 
has been required to have transverse momentum px > 200 GeV together with a pseudo-rapidity 
\rj\ < 2. For the large s region the form factor (c.f. Eq. (12)) has been used (CutOffMode = 2, 
renormScale2 = l) with t = 1. In order to simplify shape comparisons all distributions have 
been normalised to the cross section for n = 2. As expected, changing one of the two parameters 
has a small effect on the shape of the different distributions. The results are consistent, both 
in terms of shape and cross section, with previous studies. The same distributions have also 
been produced using different PDFs: MRST2001, CTEQ6L, GRV94L 7 . As expected, the effect 
of using different PDFs is very small. 

7 Although GRV94L is relatively old it has been included since it was used in some of the early graviton papers. 
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Figure 2: Differential G + jet cross section at the LHC (y/s = 14 TeV) as a function of the missing 
transverse energy. The plot shows results obtained using the form-factor compared with the case where 
all s values contribute (NoFF) and the case where the cross section is truncated (Trunc). 

Since the high energy collider signals can be sensitive to the large s region, where the effective 
theory is not necessarily valid, it is interesting to study the effect of the different methods to 
treat this region (c./. Sec. 2.3). In order to use the form factor, the quantity fi needs to be related 
to a physical scale in the production process and the choice of scale can affect the cross sections 
significantly. For the mono-jet studies presented here, the \x parameter has either been set to the 
jet energy (CutOffMode = 3) or alternatively to the jet pr (CutOffMode = 2, renormScale2 
= 1). 

Fig. 2 shows the missing transverse energy distributions from G + jet production when using 
the different options related to the UV region. The form factor results are compared with the 
case where all s values contribute (NoFF) and the case where the cross section is truncated 
(Trunc). The effect of the form factor, which implies a softer distribution, is visible together 
with the expected dependence on the model parameters. The effect decreases with increasing t, 
where t = 2 is shown to be similar to the NoFF distribution and the case with t = 0.5 approaches 
the truncated distribution. In addition, the effect increases slightly for fi = -E(jet) compared 
to ijl = pr(jet). The study showed moderate sensitivity to n and the effect becomes larger for 
smaller values of Mn, which is consistent with Eq. (12). 

3.2 Unparticle Emission 

Fig. 3 shows the differential cross section for mono-jet production due to U emission at the LHC 
as a function of the jet energy. In accordance with [7] the cross sections correspond to spin-1 
unparticle emission through the two processes, qq -^UjG + g and qg U/G + q together with 
spin-0 unparticle emission due to the gluon process, gg -^li/G + g. The distributions are shown 
for a number of dy values. The parameters = 1 TeV and A = 1 have been used for both 
the scalar and vector processes. Only jets with transverse momenta larger than 100 GeV and 
with rapidities \y(jet)\ < 2 have been considered. In addition to the cross sections predicted for 



8 




Figure 3: Differential cross section of U + jet production at the LHC (y/s — 14 TeV) as a function of 
the jet energy. Contributions from spin-0 and spin-1 unparticles are included. 

these unparticle scenarios the plot demonstrates the generally small influence of the non-integer 
dimension parameter du on the slope. Both the distribution shape and cross section values have 
been found to be consistent with the results in the reference literature. 

The coherent implementation of the similar U and G processes in Pythia8 simplifies studies 
of the possibility to distinguish between the two models. As an example, a comparison of the 
differential cross section as a function of the transverse jet energy is shown in Fig. 4. The graviton 
distributions for two values of n are compared to the distributions from spin-1 unparticle emission 
for du = 1.1, 1.5 and 1.9. The unparticle results contain both the two allowed sub-processes 
and the parameters values Ay = 1 TeV and A = 1 have been used. It is clearly visible that 
the two models exhibit different spectrum slopes almost independently of n or du and similar 
results were found for the spin-0 unparticle case. For the value, Mp = 2 TeV, the cross section 
for LED graviton production is significantly higher in the tail than for the unparticle scenarios 
shown. The cross sections, however, scale with Md or Ay which allows the cross section from 
either of the two models to be larger or smaller than the other. Therefore, as pointed out in [16], 
it is important to also study the shape of the signal in order to understand which new physics 
model is more compatible with a possible mono-jet excess at the LHC. 

4 Virtual Exchange, Photon Pair Production 

The PYTHIA8 processes related to virtual G exchange f f bar2gammagamma and gg2gammagamma 
are controlled by the program parameters LambdaT, Cut Off mode, t, and for these processes 
Cut Of fmode = 1 is not available. The corresponding unparticle processes have the same names, 
but reside in the different unparticle process category as listed in the appendix. They are 
controlled by the following parameters spinU, dU, LambdaU, lambda, with the two possible 
spin values spinU = 0,2. 
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Figure 4: Differential cross sections as a function of the transverse jet energy for LED graviton emission 
and emission of spin-1 unparticles at the LHC (y/s = 14 TeV). 

This section presents different characteristic results for di-photon production due to both 
virtual U and G exchange. The validation of these processes is based on the results presented 
in [6, 7, 9, 11, 13] and similar selection cuts as well as model parameter values have been chosen 
for easy comparison. 

4.1 Virtual Graviton Exchange 

Fig. 5 shows the invariant mass distribution of the two photons produced by LED graviton 
exchange, with and without using the form factor in Eq. (12). In accordance with [6] the mass 
scale value, Ah = 2.5 TeV, has been used 8 . The results using the form factor are shown for a 
number of different values of n. For virtual graviton exchange the renormalisation scale has been 
set to /i 2 = s (CutOff Mode = 2, renormScale2 = 4). The characteristic effects from the form 
factor are visible, in particular the softening of the invariant mass distribution. As expected 
from Eq. (12), the distributions decrease faster at large M 77 values for large n and they intersect 
at M 77 = At, which for example in Fig. 5 corresponds to 2.8 TeV. The study also showed the 
expected behaviour that the graviton contribution becomes significant at larger M 77 values with 
increased values of Ay. 

4.2 Virtual Unparticle Exchange 

The results related to 77 production at the LHC from virtual unparticle exchange have been 
compared with [9]. Cuts have been applied in all cases on the photon rapidity, | y 1 |< 2.5, 
and transverse momentum p T > 40 GeV. For the invariant mass distributions the cut on the 
polar angle, | cos(9 7 |< 0.8, has been used and for the angular and rapidity distributions only 
the invariant mass region, 600 < M 77 < 900 GeV, has been considered. Fig. 6 and 7 show 

8 The value using the Ah convention was simply translated to, At = 2.8 TeV, as described in section 2.1. 
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Figure 5: The di-photon invariant mass spectrum for LED graviton exchange at the LHC using Ah = 2.5 
TcV. The results are produced with and without using the form factor, assuming various values of n. 



the differential cross section as a function of M 77 and cos 0* respectively, for spin-0 unparticles, 
where 0* corresponds to the polar angle in the centre-of-mass frame of the photons. The 
invariant mass distribution is shown for three different A values and the signal decreases in the 
expected way with smaller coupling constants. The characteristic angular shape of the signal 
is visible and the results are found to be in good agreement with the literature listed above. 
In addition, the effects on various kinematic distributions from using different PDFs have been 
studied by using again the PDF sets MRST2001, CTEQ6L, GRV94L and the results showed 
only small differences. 



5 Virtual Exchange, Lepton Pair Production 

The available processes in Pythia8 with virtual graviton exchange leading to lepton pairs in the 
final state (Drell-Yan) correspond to f fbar211 and gg211. The relevant program parameters are 
LambdaT, CutOffmode, t, where CutOff mode = 1 is not available. Similar to the 77 case, the 
corresponding unparticle processes have the same names, but reside in the different unparticle 
process category as listed in the appendix. They are controlled by the parameters spinU, dU, 
LambdaU, lambda, gXX, gXY, with two unparticle spin values available spinU = 1,2. 

In this section we present a number of results obtained with Pythia8 used to validate the 
implementation of these processes as well as to quantify their characteristic behaviour. The 
processes are mainly validated against the results presented in [6, 7, 11, 13] and the choices 
of model parameter values and several selection criteria have also been based on these papers. 
We mainly focus the study on lepton pair production at the LHC. However, since most of the 
literature related to effects from chiral spin-1 unparticle couplings is available for e + e~ -collisions 
we also present some validation results in that context. 
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Figure 6: Invariant mass distribution of the di-photon system from spin-0 unparticle exchange at the 
LHC (y/s = 14 TeV) for different values of the A s and A t couplings. 

5.1 Virtual Graviton Exchange 

In addition to the di-photon production at the LHC, discussed in Sec. 4, the production of 
two leptons can be influenced by the effects of virtual graviton exchange. For this reason, the 
invariant mass distribution of the two final-state leptons has been studied for a large number 
of different parameter settings of the LED model and agreement with the reference literature 
has been found. As an example, Fig. 8 shows the differential cross section as a function of the 
invariant mass of the lepton pair for the cut-off scale, Ah = 2.5 TeV, in scenarios with n = 3. 
As for the 77 process the mass scale value in the common Ah convention is used. In contrast 
to the 77 case the results are given for a number of different t parameter values at a fixed 
n. The results show the strongly t-dependent increase of the cross section at high masses for 
the LED scenarios with respect to the SM. As expected from Eq. (12), the differential cross 
section converges towards the SM prediction for decreasing values of t, while for larger t values, 
e.g. t = 2.0, it approaches the LED scenarios where the form factor is not used (NoFF). In 
addition, the dependence of the cross section on the number of extra dimensions, n, has been 
studied and a similar behaviour to the di-photon production (c.f. Fig. 5) has been found as well 
as a good agreement with the results in the literature. 

5.2 Virtual Unparticle Exchange 

Similar to the virtual graviton exchange discussed above, the virtual exchange of unparticles can 
modify the invariant mass spectrum and the angular distribution of final-state lepton pairs at the 
LHC as well as in e + e~-collisions. The Pythia8 results from various settings of the unparticle 
model parameters have been compared to the original literature and again good agreement has 
been found with the results for pp collisions at the Tevatron as presented in [7] . 

Fig. 9 shows the differential Drell-Yan cross section at the LHC as a function of the invariant 
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Figure 7: Angular distribution, cos 9*, of the photons from spin-0 unparticle exchange at the LHC 
(y/a = 14 TeV). 

mass of the lepton pair for various values of du, assuming unparticles with spin-1 and vector- like 
4-fermion interactions. A cut on the lepton pair rapidity (|y| < 1) has been applied, which 
ensures that only central leptons are selected. As expected, the characteristic increase of the 
differential cross section due to the unparticle signal is more prominent for small values of du, 
e.g. du = 1.3, both above and below the Z pole. 

The possibility of chiral spin-1 unparticle couplings could imply various effects, both affecting 
the invariant mass as well as the angular distribution of the final-state lepton pair. These effects 
have mainly been studied for e + e _ -collisions in the literature and we therefore include the 
angular differential cross section obtained by Pythia8 for the SM and unparticle model, using 
various du values, in e + e~-collisions with y/s = 200 GeV. Fig. 10 shows the angular distributions 
for LL + RR (left) and LR + RL (right) interactions, where the strong impact of a chiral coupling 
on the angular distribution is demonstrated. The parameter values Ay = 1 TeV and A = 1 have 
been used. The angle 9 is defined as the emission angle of the fj,~ relative to the momentum 
vector of the incoming electron. 

The angular dependence of the differential cross section can be investigated by studying the 
forward-backward asymmetry which is commonly defined as 

N F - N B . . 

A FB = — (13) 

where Np and Nb are the numbers of forward (cos 9 > 0) and backward (cos 9 < 0) produced 
events, respectively. At hadron colliders the determination of the angle 9 is more complicated, 
since the direction of the incoming fermion is unknown. However, Apb is foreseen to be studied 
at the LHC [19] exploiting the fact that in annihilations of a valence quark and sea anti-quark the 
valence quark has, on average, a larger momentum than the sea quark. The resulting longitudinal 
momentum of the di-lepton system approximates the quark direction and, therefore, the angle 
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Figure 8: Differential cross section for the production of two charged leptons at the LHC as a function 
of their invariant mass. The results are shown for various LED scenarios with n = 3 and Ah = 2.5 TeV. 

9 between the lepton and quark in the di-lepton rest frame can be calculated. For the results 
presented here we use this method at generator level. In addition, we require, \y\ > 1, for the 
lepton rapidity and that at least one electron is in the central region, \rj\ < 2.5. 

In Fig. 11 the estimate of Afb, using the method described above, is shown for the first time 
together with effects from unparticle exchange. The two plots show the invariant mass range 
50 to 200 GeV (left) and 200 to 1600 GeV (right) respectively for various unparticle scenarios 
and the SM. It can been seen that certain unparticle scenarios (e.g. du = 1.3) lead to large 
deviation from the SM even at mass values below the Z pole which are already investigated in 
detail by past experiments [20]. However, scenarios (e.g. djj = 1.8) resulting in sizable effects 
only at larger mass values also exist, i.e. they are consistent with the SM in the already studied 
region at smaller masses. 

Experimental effects, in particular statistics, will pose a significant challenge for the recon- 
struction of Afb- A detailed experimental study is therefore needed to estimate the feasibility. 
However, in the case of an observation of a signal in the high mass tail (c.f. Fig. 8), a study 
of Afb as a function of the invariant mass could be of interest in order to confirm a non-SM 
content in the Drell-Yan spectrum. 

6 Full Event Simulation 

In order to allow for easy comparison with results in the literature as well as to clearly illustrate 
the effects related to the model parameters, the results from the various processes presented in 
this paper have been considered at the level of the hard process. 

The processes are, however, integrated within the Pythia8 framework of parton showers, 
underlying event activity, hadronisation and unstable particle decays. These steps both increase 
the complexity of the events as well as affect their kinematic properties and are normally required 
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Figure 9: Differential cross section for the production of two charged leptons as a function of their 
invariant mass. The results are shown for scenarios with unparticles using different values of du at the 
LHC. 

for proper comparisons between simulated phenomena and experimental data. 

As an example, Fig 12 shows how initial state radiation (ISR) affects the pr in G + jet 
production. At leading order the G and jet have exactly balanced pt, but additional radiation 
from ISR, or other higher order effects, breaks this balance. Fig 12 shows the pt difference 
between the G and the jet after ISR, relative to the original jet px in the hard process. This 
is shown as a function of the jet pt at the hard process level and for the different G + jet 
sub-processes separately. As seen in the plot, ISR affects the sub-processes noticeably, however, 
to different extent and the imbalance decrease for higher pt jets. This clearly illustrates the 
importance of building the processes into the simulation of the full event structure in Pythia in 
order to model such effects. 



7 Conclusions 

Both large extra dimensions as well as unparticle models have become very popular among pos- 
sible extensions of the SM. Many of the main processes related to LED searches have analogous 
processes in the unparticle scenario, which are often more general in terms of model parameters. 
A common way to implement the similar processes has been found within the Monte Carlo 
generator PYTHIA8, where only a few translations of model dependent constants are necessary. 
In addition to technical benefits, the common implementation makes comparisons between the 
processes from the two models consistent and simplifies studies of how differences at the process 
level for a certain final state can affect the measured observables. The processes can be used 
together with several options for the treatment of the UV region of the effective theory. In this 
paper, the options related to the individual processes have been discussed together with results 
obtained by Pythia8. In addition to providing references for the Pythia8 program, the results 
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Figure 10: Angular distribution for e + e — > including unparticlc exchange, for various values of 

du at y/s = 200 GeV, assuming LL + RR (left) and LR + RL (right) 4-fermion interactions. 
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Figure 11: as a function of the invariant mass of the lepton pair at the LHC. The results are shown 

for du — 1.3, 1.8 with vector-like and LL + RR unparticle interactions (left) as well as for du = 1.8 with 
various chiral interactions (right). In both figures spin-1 unparticle exchange is assumed. 



illustrate characteristic effects of the model parameters and various key features of the different 
kinematic distributions. The results presented in this paper have been focused on mono-jet, 
di-photon and di- lepton final states at the LHC. 
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Figure 12: Relative pr difference between the graviton and jet due to initial state radiation. The results 
are shown for G + jet production with n = 6 and Mp = 4 TeV. 

Appendix -U/G Specific Parameters in Pythia8 

In the Pythia8 program the U/G processes are turned on/off by the following parameters, 



ExtraDimensionsUnpart : f f bar2Ugamma - Real 14 emission, // — >■ 14 + 7. 

ExtraDimensionsUnpart : f f bar2UZ - Real 14 emission, // — > 14 + Z . 

ExtraDimensionsUnpart :qqbar2Ug - Real IA emission, qq — > U + g. 

ExtraDimensionsUnpart : qg2Uq - Real IA emission, qg — > U + q. 

ExtraDimensionsUnpart : gg2Ug - Real IA emission, gg — ^ 14 + g. 

ExtraDimensionsLED : f fbar2Ggamma - Real G emission, // — >■ G + 7. 

ExtraDimensionsLED : f fbar2GZ - Real G emission, f f G + Z. 

ExtraDimensionsLED : qqbar2Gg - Real G emission, qq — > G + g. 

ExtraDimensionsLED : qg2Gq - Real G emission, qg — > G + q. 

ExtraDimensionsLED : gg2Gg - Real G emission, gg — >■ G + g. 
ExtraDimensionsUnpart : ffbar2gammagamma - Virtual IA exchange, // — > (IA*) 

ExtraDimensionsUnpart : gg2gammagamma - Virtual IA exchange, gg — > (14*) - 

ExtraDimensionsUnpart : ffbar211 - Virtual IA exchange, // — > (14*) 

ExtraDimensionsUnpart :gg211 - Virtual IA exchange, gg — > (14*) - 

ExtraDimensionsLED : ffbar2gammagamma - Virtual G exchange, // — > (G*) 

ExtraDimensionsLED : gg2gammagamma - Virtual G exchange, gg —> (G*) - 

ExtraDimensionsLED :ffbar211 - Virtual G exchange, // -)• (G*) 

ExtraDimensionsLED :gg211 - Virtual G exchange, gg — > (G*) - 



■ 77- 
77. 

■II. 
It. 

> 77. 
77. 

> it. 

el 



The model parameters related to the LED graviton processes are specified by, 
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ExtraDimensionsLED:n - Number of large extra dimensions, n. 

ExtraDimensionsLED:MD - Fundamental scale of D dimensional gravity, Mp. 

ExtraDimensionsLED:LambdaT - Cut-off scale for virtual G exchange, Ay. 
ExtraDimensionsLED : CutDff mode - This parameter specifies the treatment of the high 

energy contributions. Possible values are to 3. 
ExtraDimensionsLED :t - Form factor parameter, t. 



The model parameters related to the unparticle processes are specified by, 



ExtraDimensionsUnpart : spinU 
ExtraDimensionsUnpart : dU 
ExtraDimensionsUnpart :LambdaU 
ExtraDimensionsUnpart : lambda 
ExtraDimensionsUnpart :gXX 



ExtraDimensionsUnpart :gXY 



ExtraDimensionsUnpart : CutOf f mode 



The unparticle spin. Possible values are, 0, 1 or 2. 
Scale dimension parameter, djj. 
Renormalization scale, Ay. 
Coupling to the SM fields, A. 

Options for the LL/RR helicity contributions, gxx- 
This parameter is only relevant for the spin-1 11 
process, where the available options 0, 1 and 2 
correspond to the helicity amplitudes being 
multiplied by a factor 1, or -1. 
Options for the LR/RL helicity contributions, gxY- 
This parameter is only relevant for the spin-1 U 
process, where the available options 0, 1 and 2 
correspond to the helicity amplitudes being 
multiplied by a factor 1, or -1. 
This parameter specifies the treatment of the high 
energy contributions. Possible values are or 1. 



More details on the meaning of the model parameters are given in the text. 

Appendix - Scalar Unparticle Production, // — > ZU. 

The relevant scalar unparticle interaction terms for this process are given by, 

£ D ^^/(A 1 + A^75)/ (14) 



A 



Due to CP conservation, the two couplings Ao or A' are only considered separately. Since the 
matrix element expression of the two cases becomes identical in the absence of any interference 
with the SM, only one common A value is used for the implementation in accordance with Eq. 
(3). These couplings are usually suppressed by a quark mass factor, because of the chirality 
flipping, and are therefore expected to be small in the scalar unparticle case. However, in the 
phenomenological spirit of this paper the choice of restricting model parameter values is left to 
the user. 

As described in [5], the cross section for the, // — > ZU, process is given by, 



d 2 a |M| 2 A du (dP l 



dPfidt 16vrs 2 2ttA 2 I A 2 



2 \ d-u-2 

^' (15) 
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where \M\ is the colour and spin averaged matrix element and is the invariant mass of the 
unparticle, < Pfi < (^/s — mz) 2 - The matrix element is given by [21], 



{gt q + g% q )AA\ 2 




(16) 
(17) 
(18) 



with g Ljd = -\ + \sm 2 9 w , g R4 = +±sm 2 6 w , 9l,u = \~ § sin 2 9 W and g R>u = -fsin 2 ^. 
Regarding the amplitude formula, \A\ 2 , the first two terms inside the square brackets result 
from the t-channel, the following two terms from the u-channel and the last term from the 
interference. 

In the same way as for the other spin versions of this process, the process where a photon is 
emitted together with an unparticle has been obtained by the following photon limit of the Z 
matrix element, 



9l±9r 



mz —s- 

>Q 2 
-> l 



sin 6w cos 9w 

where Q is the electric charge of the incoming fermions. 



(19) 
(20) 

(21) 
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